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Low energy electron diffraction (LEED) experiments, LEED simulations and finite slab density 
functional calculations are combined to study the cleavage surface of Co doped BaFe2-a;CoccAs2 (x 
= 0.1, 0.17). We demonstrate that the energy dependence of the LEED data can only be understood 
from a terminating 1/2 Ba layer accompanied by distortions of the underlying As-Fe2-As block. As 
a result, surface related Fe states are present in the electronic structure, which we identify in 
angle resolved photoemission experiments. The close proximity of the surface-related states to the 
bulk bands inevitably leads to broadening of the ARPES signals, which excludes the use of the 
BaFe2-ccCoccAs2 system for accurate determination of self-energies using ARPES. 



The recently discovered iron pnictide superconductors 
[1] have quickly turned into one of the most widely stud- 
ied systems in condensed matter research, using all the 
modern tools that have been perfected through years of 
research on another class of high temperature supercon- 
ductors: the cuprates. However, two of the most impor- 
tant and widely used tools to determine the electronic 
structure, scanning tunneling microscopy (STM) and an- 
gle resolved photoemission (ARPES) are surface sensitive 
techniques. It is therefore crucial to assess the impact of 
the surface on the electronic structure. 

Several different structures containing As-Fe2-As 
blocks are known: LnOFeAs, MFeAs and MFe2As2 com- 
monly referred to as the "1111", "m" and "122" fami- 
lies. Subdividing the unit cell into As-Fe2-As blocks and 
spacer layers results in an alternating stack of positively 
and negatively charged layers. Because of the alternat- 
ing charges of these layers, the presence of a surface can 
give rise to a diverging electric potential in the solid 
[2, 3 . Such polar surfaces are energetically unfavorable 
and usually a structural or electronic reconstruction takes 
place to prevent their occurrence. LaOFeAs is such an 
example: cleavage takes place between the LaO and As- 
Fe2-As blocks, resulting in either a highly reconstructed 
La or As termination [4 which supports surface-related 
states that have been observed in ARPES [5]. The sur- 
face of the 111 system, on the other hand, barely influ- 
ences the electronic structure [6^. 

This leaves the MFe2As2 pnictide superconductors 
which are by far the most widely studied using surface 
sensitive probes. In this case, STM experiments show a 
diversity of topographies resulting from the cleaving pro- 
cess [71412]. which have been interpreted as either con- 
sisting of As or Ba layers. Recent first principle calcu- 
lations for the parent compounds show a reconstructed, 
Ba terminated surface to be energetically favorable, leav- 
ing open the possibility for a metastable As termination 



upon low-T cleavage ^3]. There is a paucity of well- 
founded experimental data on the exact surface struc- 
ture of the superconducting 122 compounds, in particu- 
lar when cleaved at cryogenic temperatures as is usual 
for ARPES and STS studies. Consequently, the impact 
of the real surface structure on the electronic states re- 
sponsible for superconductivity remains an open and im- 
portant issue. 

In this Letter we combine IV-LEED experiments and 
simulations with DFT slab calculations to investigate the 
cleavage surface of superconducting, Co-doped Bal22. 
The data show that the half Ba-layer termination has 
an impact on the underlying As-Fe2-As block, leading to 
significant departures from the bulk structure. This gives 
rise to surface-related Fe Zd bands, which, in turn, can be 
found in our ARPES data. These surface states impede 
accurate quantitative analysis of (for example the widths 
of) ARPES features in the Bal22 systems. 

We start with the LEED experiments. Samples were 
top-post cleaved in-situ at pressures better than 5xl0~^^ 
mbar at temperatures under 20K - following the standard 
recipes used in ARPES and/or STS studies. For details 
on samples and experimental procedures, see supplemen- 
tary materials. The measured LEED patterns for x = 0.1 
and 0.17 samples [Fig. [TJa),[TJb)] resemble a 2 x 2 pat- 
tern, which in fact arises from a combination of ordered 
surface structures as indicated in Fig. [ijc). In order to 
circumvent the polar surface (type III in the classifica- 
tion of Ref. [2 ), half of the Ba layer goes with each half 
of the cleaved crystal. This now places a terminating 1/2 
Ba layer with charge +1 on top of the underlying As-Fe2- 
As block which has overall charge -2. We obtain a zero 
net dipole moment if we think of the Ba layer below the 
uppermost As-Fe2-As block as half belonging to this up- 
per block and half belonging to the next As-Fe2-As block 
below. Continuing in this fashion represents the only en- 
ergetically favorable situation of any likelihood [TT. Or- 
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FIG. 1: (Color online) (a): Typical LEED images for Eo = 
114 eV and 292 eV for BaFe2-xCo^As2 with x = 0.17. (b) 
Eo = 200 eV for X = 0.1. In each image, representative (1,1) 
spots are highlighted with red circles. All images taken at 
17 K. (c): real space and q-space structures of the 2x1 and 
y^2) X \/(2) structures (green: Ba; purple: As; orange: Fe) 
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FIG. 2: (Color online) Half integer IV-curves for both termi- 
nations. Left: 2x1 termination, right: V2 x a/2 termination. 
The Pendry R factor is indicated for each spot. 



dering in the partial Ba termination layer gives rise to the 
additional, fractional spots in the LEED images as indi- 
cated in Fig. [ijc). The patterns observed are identified 
as 2 x 1, 1x2 and a/2 x a/2 ordered partial Ba termi- 
nation layers. The integer 1x1 pattern contains higher 
order reflections from all types of superstructure. The 
reconstructions are not limited to the tetragonal phase 
only: the LEED pattern of the underdoped, orthorhom- 
bic phase is identical, as shown in Fig. [ijb). 

Even though the arguments presented above on the ba- 
sis of energetics are already persuasive, we now prove the 
1/2 Ba termination using IV-LEED data, in which we 
solve for the real surface structure of the outermost four 
atomic layers. Figure [2] presents the results from an x = 
0.17 sample, together with the relevant simulations, car- 
ried out using the Barbieri/Van Hove SATLEED package 
[15] (see supplementary material for details). Pendry R- 
factors [16] are used to measure the level of agreement be- 
tween experiment and simulation. Simulations assuming 
a bulk crystal terminated by a 1/2 Ba layer with either 
\/2 X a/2 or 2 X 1 superstructure resulted in overall Pendry 
R-factors, Rp ~ 0.31 and Rp ~ 0.46, respectively. As a 
satisfactory solution should have Rp < 0.3, the ordered 
superstructures of the half-Ba layer alone are not suffi- 
cient. We therefore allowed the atomic positions of the 
deeper lying layers to relax. Given the measured energy 
ranges (in total 1208 eV for the x and 1782 eV for 
the 2x1 superstructure), optimizing the parameters of 
the first four atomic layers so as to give reliable estimates 
of the actual deviations from bulk positions proved to be 
feasible. We arrive at final Pendry R-factors of Rp = 0.19 
for the resulting distorted a/2 x a/2 termination and Rp 
= 0.29 for the distorted 2x1 structure. IV-simulation 
for a 2 X 2 reconstruction was unable to arrive at Rp < 
0.5 for any 2x2 structures we adopted. We also tried As 



terminated surfaces, also relaxing the atomic positions 
in the first three layers. This gives much higher Pendry 
R-factors, Rp ^ 0.48 for the 2x1 structure and Rp ^ 
0.42 for the a/2 x a/2 structure. Therefore, based both 
on robust electrostatic considerations and our quantita- 
tive IV LEED study we can conclusively state that the 
termination surface of the Co-doped Bal22 superconduc- 
tors cleaved at cryogenic temperatures consists of a half 
Ba layer coupled to relaxed atomic positions in the near- 
surface layers. 

The atomic positions from IV-LEED are given in ta- 
bles in the supplementary information, and clearly indi- 
cate significant deviations from the bulk structure. Apart 
from a small reduction in the inter-planar distances, we 
also observe buckling of the outermost As-Fe2-As block 
and lateral displacements of the atoms in the first four 
layers. The Debye temperature of the surface Ba atoms 
is found to be quite low (60 K), making these Ba atoms 
relatively mobile. This matches the observation that the 
fractional spots irreversibly fade out when the sample is 
warmed to approximately 150 K [10]. STM experiments 
also report destruction of the long range order seen in 
the topographs by thermal cycling to high temperature 
[7 . Incidentally, this also explains why no clear frac- 
tional spots where observed in Ref. [9 , as those samples 
were cleaved at room temperature. Seeing as almost all 
ARPES studies are conducted on cryogenically-cleaved 
samples, the low T structure reported here is of direct 
relevance. 

Now that we have proven the existence of significant 
(sub) surface reconstructions, the next step is to quan- 
tify their effects on the electronic structure. This we do 
via density functional theory (DFT) calculations for the 
different surface geometries using the full potential local 
orbital code (FPLO) [17^ (for details see supplementary 
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FIG. 3: (Color online) Top view (a) and side view of half a 
slab (b) used in the DFT calculations (green: Ba; purple: As; 
orange: Fe). The blue shaded area in (b) indicates distorted 
atomic positions, as obtained from the LEED experiments. 



information) . We simulated the cleavage surface by using 
z-periodic replicas of a finite slab, separated by vacuum 
(i.e. repeated slab calculations). The structure of the 
slabs used in the calculations is shown in Fig. [sj^a). In 
Fig. [sj^a) and (b) we only show the top half of the slab: 
the other half is identical, but mirrored in the lowest Ba 
layer shown. To distinguish the effect of the Ba super- 
structure from the distortions, we repeated the calcula- 
tions with similar slabs that have all atoms at the bulk 
positions, but with a 2 x 1 or >/2 x >/2 ordered half-layer of 
Ba at the surface. To address doping, the virtual crystal 
approximation (VGA) is employed. The Ba superstruc- 
tures impose a planar super cell on the whole slab, which 
has the consequence that all bands, including the bulk 
bands, are folded back into the smaller Brillouin zone 
of the super cell. Obviously, this back folding becomes 
less and less meaningful as the layers become increasingly 
bulk-like. In the following, we unfold the band structures 
according to the translational symmetry of the original 
planar unit cell using the approach of Ref. [T8^ (see sup- 
plementary information for details). Figure [4]^ a) shows 
an experimental FS for the x = 0.17 sample measured 
at 74 eV photon energy. The DFT bulk FS is shown 
in Fig. ^h). Comparing the two, we note two distin- 
guishing features: (i) a weak FS sheet at the X point 
in the experimental data and (ii) the fact that the elec- 
tron FS around M breaks the fourfold symmetry of the 
bulk crystal. The first feature has also been observed 
in SrFe2As2 [7], and disappeared when the sample was 



FIG. 4: (Golor online) (a): Fermi surface of 
BaFei.83Goo.i7As2 measured with 74 eV, p-polarized 
photons. Indicated in yellow is the 1st Brillouin zone and 
the r and M points. The dashed yellow line indicates the 
cuts shown in Fig. [s] (b): undoped bulk FS from DFT. 
(c): Unfolded FS for the 2x1 termination, from DFT using 
the relaxed atomic positions and VGA doping of 0.17. Red: 
bulk; green: surface; brown: mixture of bulk and surface 
states, (d): same as (c), but for the a/2 x y/2 structure. 



heated to 200 K. The same feature is also in the Gal 22 
data of Ref. [19,, making it a ubiquitous feature of the 
M122 iron pnictide family. Comparing with Fig. [4]^c), 
we see that this feature arises from the 2x1 and 1x2 
superstructures. At this stage, we cannot decide which 
out of the 2 X 1 or \/2 x \/2 structures is responsible for 
(ii), as both reconstructions lead to significant deviations 
from four-fold symmetry. We also note that we do not 
observe the large distortion-induced FS pocket encircling 
the r point predicted by the calculations. More extensive 
variation of the experimental ARPES parameters (pho- 
ton energy, polarization, geometry) will be necessary to 
judge to what extent this is a matrix element effect. 

In Fig. |5) we compare an /(/c, £^) image taken along the 
T — M direction in panel (a) with the bands calculated for 
the (c) 2 X 1 and (d) ^/2x^/2 distorted structures. As has 
been observed in all comparisons between ARPES and 
DFT calculations in the iron pnictides, one can identify 
most features if one includes a scaling factor renormaliz- 
ing the energy scale of the DFT data by a factor of 2-3. 
Apart from this re-scaling, there is good agreement be- 
tween the calculated and measured spectra. In particular 
we observe a surface state, indicated in Fig. [sja) by the 
dashed yellow curve that gives rise to the central peak 
in the displayed EDG. Further proof for the identifica- 
tion of this band with surface-related states comes from 
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FIG. 5: (Color online) (a): /(/c, E) image taken along the cut 
indicated in Fig. |4ja). On the side we show an EDC taken 
along the line indicated by the arrow. The dashed curve in- 
dicates the surface state (b): same as (a) but taken after 
cycling in temperature, (c) and (d): calculated band struc- 
ture for the same momentum window. Note that the energy 
window is twice as large for the DFT data. 



Fig. jsjb), which is the same cut also measured at 10 K, 
but taken after cycling the sample temperature to 150 
K. Like the fractional spots in the LEED experiments, 
the surface-related features fade away on the irreversible 
melting of the ordered Ba superstructure. Comparing the 
complexity of surface and bulk bands in panels (sjc) and 
[sjd) with the experimental data, several surface states 
probably contribute to the experimentally observed sur- 
face state. The closeness of the surface and bulk states 
will inevitably lead to a broadening of the ARPES spec- 
tra. It is likely that most of our conclusions carry over 
to other members of the 122 family, as both termination 
surfaces have been observed with STM. Cal22 appears 
to be somewhat special due to the fact that there is a 
larger energy separation between the 2x1 and a/2 x a/2 
termination states [13 . This is also evident from STM 
experiments where only 2x1 termination is observed [20 . 

In conclusion, our combined IV-LEED (experiment 
and theory) and band mapping (experiment and the- 
ory) have clearly shown the surface of the Bal22-based 
high temperature pnictide superconductors to be not 
only comprised of a 1/2-Ba layer, but also to be sig- 
nificantly distorted with respect to the bulk structure. 
This gives rise to surface-related bands, which can be 
clearly identified in ARPES data. These conclusions go 
a long way to resolving the riddle as to why so much 
of the published ARPES data on the Bal22 family is 
so broad: this is due to closely overlapping surface and 
bulk-derived states. This essentially excludes the use of 
the Bal22 system for accurate determination of in par- 
ticular the imaginary part of the self-energies in these 
superconductors via ARPES. 
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SUPPLEMENTARY MATERIAL 



EXPERIMENTS 

Crystals used in the experiments where grown at the 
University of Amsterdam from self flux. All samples are 
characterized by resistivity measurements. The critical 
temperatures are Tc = 14 K and 24 K for x=0.10 and 0.17 
samples respectively. The x=0.10 sample has a struc- 
tural/magnetic transition around T^^ = 43 K and T^ag 
= 35 K. 

The LEED experiments where performed using a mod- 
ified ErLEEDlOOO setup from Specs that allows for si- 
multaneous beam-current readout. The spot size for this 
setup is < 1 mm. The cleaved crystals were aligned with 
the surface normal axial to the LEED set-up to within 
0.2 degrees. The images where captured using a stan- 
dard CCD camera at 2 eV intervals in primary beam 
energies. For each set of experiments a dark field image 
was subtracted from the data. To account for changes 
in background intensities, we subtracted the background 
intensity averaged over a square box around the main 
spot as detailed in [21]. 

The ARPES experiments where performed at the SIS- 
HRES beamline of the Swiss Light Source, using a Sci- 
enta R4000 analyzer with a total energy resolution (pho- 
ton + electron) of 16 meV (at T = 10 K, hz/ = 74 eV) 
and a momentum resolution better than 0.05 7r/a along 
the slit. 



LEED calculations 

We have performed two distinct LEED calculations, as 
the two coexisting structures (2x1 and >/2 x ^/2) gen- 
erate different beam sets. Only half-integer beams have 
been used since the integer beams are created from the 
intensities of several structures at once. The fractional 
order beams belong to one of the calculated structures 
provided that there are no ordered structures with cell 



TABLE I: The exact coordinates in the a/2 x a/2 termination. 
The bulk values (z,x,y) are quoted in brackets. 



Ba 


0.1717 


-0.1834 


-0.2012 


0.0000 


0.0000 


0.0000 




2.0011 


1.5394 


1.6523 


As 


(1.8926 
2.0261 


1.9780 
1.5059 


1.9780) 
-2.1909 




(1.8926 


1.9780 


-1.9780) 




3.2869 


3.9686 


-2.1561 




(3.2363 


3.9559 


-1.9780) 




3.3294 


3.7308 


1.7104 


Fe 


(3.2363 
3.3778 


3.9559 
5.7287 


1.9780) 
-0.3662 




(3.2363 


5.9339 


0.0000) 




3.3832 


1.8588 


-0.3138 




(3.2363 


1.9780 


0.0000) 



size greater than 2 times the bulk unit cell. We have as- 
sumed that there are domains in every possible orienta- 
tion, i.e. four different orientations for both surface struc- 
tures, and also that there are equal quantities of 2 x 1 and 
1x2 domains. Since all orientations are observed with 
roughly equal frequency in STM experiments on a single 
cleave, we average the experimental IV curves of sym- 
metry equivalent spots (eg. (1,1), (1,-1), (-1,1), (-1,-1) are 
averaged). In the calculations no symmetry is enforced 
and equivalent spots have been averaged. As a result 
the calculated LEED patterns have four fold symmetry 
even if the surface structures do not have this symmetry. 
The bulk unit cell of this sample is tetragonal so that the 
surface unit cells are exactly 2x1 and ^/2 x ^/2. Resis- 
tivity measurements on the crystals used show no sign of 
the tetragonal to orthohombic transition for the x = 0.17 
sample for which the calculations where performed. The 
tetragonal two dimensional unit cell has been assumed 
to be of size 3.9624 A x 3.9624 A. This is an average of 
quoted values [22l |23] . The LEED calculations were per- 
formed using the Barbieri/ Van Hove SATLEED package 
[15] . The 12 phase shifts needed in the LEED calcula- 
tion were calculated with the Barbieri/ Van Hove phase 
shift package [24]. Temperature effects were calculated 
within the SATLEED code by multiplying each atom's 
scattering amplitude by a Debye- Waller factor. Pendry 
R-factors [16] were used to measure the level of agreement 
between measured and calculated IV spectra. Statistical 
errors in analysis were estimated with Pendry RR-factors 
[T6] . We have used the averaged T- matrix approxima- 
tion (ATA) |25J to take into account the Co-doping. In 
this approximation the Fe/Co -atoms are replaced by an 
averaged atom, with a scattering amplitude which is a 
weighted average of the scattering amplitudes of Fe and 
Co. The results would not differ much if only Fe would be 
used in this place as the Fe and Co are very similar scat- 
terers and the amount of Co is not too big. During the 
calculations the surface geometry and the complex inner 
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TABLE II: The exact coordinates in the 2x1 termination. 
The bulk values (z,x,y) are quoted in brackets. 



SDdb 


0.2730 


-0.0243 


-0.3503 


0.0000 


0.0000 


0.0000 




2.0895 


1.8537 


1.5570 


As 


(1.8926 
2.1125 


1.9780 
5.7112 


1.9780) 
1.7153 




(1.8926 


5.9339 


1.9780) 




3.3173 


1.9307 


-0.1066 




(3.2363 
3.3266 


1.9780 
6.0300 


0.0000) 
-0.3240 


Fe 


(3.2363 
3.3885 


5.9339 
3.7040 


0.0000) 
1.9820 




(3.2363 
3.4274 


3.9559 
-0.1746 


1.9780) 
2.1709 




(3.2363 


0.0000 


1.9780) 



potential were optimized. The Debye temperatures have 
been optimized for each element. The Ba surface- atoms 
have a different Debye temperature from the ones deeper 
in the bulk. The Debye temperature for Fe is optimized 
for Fe solely. Because the of Co and Fe are close to 
each other and the amount of Co is small, the Co doping 
should not affect the optimized Tjj for Fe. 

The optimized atomic positions obtained for the two 
different reconstructions are collected in table [J for the 
^/2 X terminated surface and in table [ll] for the 2x1 
termination. Compared to the bulk structure [22 , sub- 
stantial lateral changes occur in the three topmost lay- 
. The Fe-layer approximately 3 A below the 



III 



ers [table 

surface Ba-atoms still has lateral changes of the order 0.4 
A in both terminations. However, since LEED is more 
sensitive to geometric parameters normal to the surface 
than lateral parameters, the error margins on the quoted 
lateral movements are quite big. The estimated errors in 
lateral geometric results are given for the direction join- 
ing the optimized atomic position and the bulk position. 
This means that the Pendry RR-factors are calculated 
along the line from bulk position to the best fit position. 
The errors are calculated one at a time, while the other 
atoms reside on the best fit positions. On doing so, all 
lateral coordinates have error margins of the order ±0.2 
A and all coordinates normal to the surface have error 
margins of the order ± 0.04 A. Apart from the lateral 
movements the first As-Fe2-As block below the terminat- 
ing Ba layer becomes buckled and the interlayer distances 
dimisnish. 

When integer beams and half integer beams belong- 
ing to the a/2 X a/2 termination are taken simultaneously 
into the calculation, one gets a structure that is nearer 
to the bulk structure but the total Pendry R-factor in- 
creases to 0.32 from the 'only half integer beams' -value 
0.19. The dz- values grow towards the bulk values (1.84 
A,1.32 A, 1.27 A) and buckling diminishes (O.OlA for the 
As layer, 0.07A for the Fe layer). 



TABLE III: The geometric structure of the BaFe2-ccCoccAs2 
surface for the two different surface structures. The distances 
dzBa-As, dz As-Fe and dzFe-As are measured between the 
highest and lowest atom in neighboring sheets. Average error 
bars on the inter-planar distances are dz 0.04, while those on 
the lateral movements are d= 0.18. All values are in A. 



bulk V2xV2 2x1 



dZBa-As 


1.89 


1.83 


1.82 


dZAs-Fe 


1.34 


1.26 


1.20 


dZFe-As 


1.34 


1.22 


1.20 


buckling 








As- layer 




0.03 


0.02 


Fe-layer 




0.10 


0.11 


lateral shifts 








B a- layer 




0.27 


0.35 


As- layer 




0.55,0.52 0.44,0.34 


Fe-layer 




0.18,0.35 0.12,0.34 



0.42,0.34 0.26,0.26 



The thermal movement of the atoms is taken into cal- 
culation by fitting Debye temperatures. For 2x1 termi- 
nation, the best fit value for Ba is 60i^ which is quite low 
compared to the tabulated bulk value of llOi^T [26]. For 
a/2 X a/2 termination the best fit value is 70K. These val- 
ues are lower than the usual bulk value divided by a/2 for 
surface atoms which would give 78i^. It means that the 
surface Ba is very mobile which is to be expected since 
it can form different terminations at low temperatures. 
For As the tabulated value for the Debye temperature 
is 282i^. We find that the Debye temperature is much 
higher than that but the sensitivity in the simulations to 
the Td{As) is so low that no precise value can be given. 
The Td for Fe is 675K for V2x a/2 termination and 425K 
for 2 X 1 termination. These values are of the same order 
of magnitude as the tabulated value for Fe (470i^). 



LDA calculations 

We use the experimental lattice parameters and Wyck- 
off positions and focus on temperatures above the spin 
density wave transition, which allows for non-spin polar- 
ized calculations. Given, that we do not relax the lattice 
we also do not have to worry too much about the choice 
of the exchange and correlation functional used in the 
calculations. Experience shows, that although lattice de- 
grees of freedom are notably influenced by the choice of 
the exchange and correlation functional, band structures 
are quite insensitive to the choice (given the same lattice 
structure). 

In order to correctly treat the boundary condition of 
the Poisson equation we used a symmetric slab, with a 
cr^-mirror plane in the center of the slab. In that way we 
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FIG. 6: The bulk bandstmcture of undoped BaFe2As2. The 
band structure is colored according to orbital weights. 
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FIG. 7: Demonstration of the band unfolding technique for 
the doped ordered 2x1 slab, a) Band structure without 
unfolding, b) unfolded band structure. c,d) Fermi surface 
without and with unfolding. 



prohibit the building up of artificial surface dipoles. The 
self consistent calculations used a 6 x 6 x 1 /c-mesh and the 
tetrahedron method for the Brillouin zone integration. 
Note, that the slabs have larger in-plane unit cells than 
the bulk setup, due to the surface Ba superstructures. 
That makes the number of /c-points taken sufficient. Rel- 
ativistic effects were treated on a scalar relativistic level. 
We use repeated slab calculations in which the slabs are 
separated by vacuum. We made sure that the vacuum 
between the slab is wide enough to decouple the repli- 
cas. In a slab calculation the /c^-dispersion is basically 
replaced by a projection of the band structure onto the 
kx^ky-'^Xd.ne. Due to the restriction of the calculations 
to a rather small finite slab we will only get a limited 
number of Fermi surfaces in the region of the projected 
A:;2-dispersion. In our case we have two bulk and two sur- 
face FeAs triple layers, which is of course not enough to 
show all of the /c^ -projected bands. However, we will get 
representatives of the bulk bands in the correct regions in 
/c-space, which can be compared to pure surface bands. 

There are several indicators for convergence of a finite 
slab to a finite slab with bulk like electronic structure 
far away from the surface. The first is, that the inner 
most layers of the finite slab have a population analysis 
(and hence band structure/density) unchanged from the 
bulk calculation. The slabs we used have four Fe2As2 
triple layers and it turns out that the two inner triple 
layers agree in their population numbers with the bulk 
results within a few hundredth of electrons. The second 
issue is that a charged surface will locally change the 
chemical potential, which manifests itself in band shifts 
of the surface bands compared to the bulk bands. This 
difference of the Fermi level between the surface and the 
bulk will gradually vanish going from the surface into the 
bulk. Hence, another parameter to be converged is the 
Fermi level of the bands belonging to the inner most lay- 



ers. In our case the surface charging is less pronounced 
than in the 1111 case [4 since the formation of Ba su- 
perstructures restores the local stoichiometric chemical 
composition, which in turn makes our four triple layer 
slab super cells large enough. 

To illustrate the change in occupations, we focus on 
the Mullikan type population analysis obtained by pro- 
jection onto the local orbitals. Due to hybridization and 
overlap effects the occupation numbers do not reflect one 
to one the chemical valency, however trends can be ob- 
tained quite easily. In the bulk the Ba misses one elec- 
tron, while the two As per formula unit each have 0.55 
excess electrons (while the difference to neutrality sits on 
the two Fe atoms). In the slab calculations the single 
surface Ba misses only 0.49 electrons, while the top most 
(sub surface) As layer now has 0.33 excess electrons on 
average. In the surface superstructures considered here 
every second surface Ba atom is missing while the FeAs 
triple layers have two Fe2As2 formula units. Counting 
Ba atoms correctly reveals that each four-Fe triple layer 
must be combined with one Ba from above the layer and 
from below the layer in order to obtain the stoichiomet- 
ric BaFe2As2 formula unit. So the two formula units per 
unit cell of the surface layers have 1.5 missing electrons, 
which when compensated by the top most As layer leads 
to an As occupation |, which is quite close to what we 
obtained. The situation is more or less the same for all 
slabs considered, although there are some variations in 
how much the charge compensation sits only on the As 
layer under the surface compared to also sitting to a lesser 
degree on the As layer below the topmost Fe layer. 

In a recent publication [18] a technique for band un- 
folding was described. The original idea is based on Wan- 
nier functions, however as is pointed by the authors a 
local orbital basis is a suitable starting point as well. In 
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FIG. 8: (Color Online) Comparison of slabs with superstructure only ((a),(c)) and with distortions in the Ba layer and As-Fe2-As 
block ((b), (d)). The top (bottom) panels are for the 2x1 (v^ x a/2) superstructure. 



fact it turns out that the technique is even simpler in 
a local basis code. The usual "fat-band" plots are ob- 
tained by projecting the Bloch states onto local orbitals 
of a certain type and symmetry, for instance Fe 3d or- 
bitals. One can spin the symmetry selective properties 
of such a projection a bit further by including transla- 
tional symmetry, which is achieved by using Bloch sums 
of these local orbitals as projectors. If the phase fac- 
tors of the Bloch sums are chosen such that they become 
representations of a certain translational symmetry, the 
resulting band weights will show destructive interference 
for states, which are replicas of the bands according to the 
chosen translational symmetry. Ref [18] points out that 
the Bloch-sum orbital weights as obtained by this ap- 
proach contain some of the matrix element effects, which 
enter real life measurements such as ARPES. We imple- 
mented this technique into FPLO and will use it in what 
follows in order to simplify the discussion. The trick al- 
lows not only the unfolding of bands but naturally also 
includes projection on certain layers of the slab. 

Figure |6] shows the well known band structure, while 
the Fermi surface of undoped bulk BaFe2As2 for T > T^t 
> T^mag IS showu in Fig. 4(b) of the main paper. The 



line thickness in Fig. [6] represents the orbital character 
of the corresponding bands. The main contributions to 
the bands around the Fermi level comes from the Fe 3d 
orbitals (red), with small admixtures of As 4p (green). 
There is a special feature of the Bal22 compound be- 
tween and 400 meV above the Fermi level along the 
line FX, which is due to a Ba 5d band. Since, the iron 
character is so dominating we will use the Fe 3d orbitals 
in the following, whenever projections onto surface and 
bulk layers are needed. 

We have performed calculations for four different slab 
geometries, i.e. the two distorted slabs as shown in Fig. 
3 of the main paper and the corresponding undistorted 
(ordered) slabs, which have the a/2 x a/2 and 2x1 pla- 
nar Ba surface superstructures but bulk positions and 
distances otherwise. 

The consequences of the unfolding technique are illus- 
trated in Fig. [7| which shows the band structure (BS) 
and Fermi surfaces (FS) for the undistorted 2x1 slab 
along the symmetry lines of the simple planar BZ of the 
Fe2As2 unit cell. In the upper panel the standard band- 
structure calculation result without unfolding is shown. 
The backfolding due to the 2x1 super cell is clearly seen 
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FIG. 9: Unfolded Fermi surfaces of the 2 x 1 (left) and v^x 
(right) Ba surface layer super structures. Upper (lower) panel: 
with undistorted (distorted) surface and sub-surface layers. 
The unfolding technique was used. Green: bands from the 
Fe-surface layer. Red (brown): bulk (mixed) states. 



along the line FX. This is also illustrated by the FS in 
Fig. [tJ^c), where replicas of the FS sheets around F are 
found around the X-point. The green sheets are predomi- 
nantly of Fe surface layer character (the iron layer closest 
to the surface), while the red (brown) sheets are predom- 
inantly of bulk (mixed) character. It is easy to under- 
stand that the backfolded bands around X have nearly 
exclusively surface character, since the effect of the Ba 
surface superstructure on the electronic structure decays 
with increasing layer depth. A similar discussion holds 
for the replica of the bands around the M-point, which 
appear at Y(M'). In Fig. ^h) and [7]^d) the results af- 
ter unfolding are shown. We chose the same color code. 
Most notably, the back folded bands around the X-point 
fade away due to the destructive interference resulting 
from the unfolding projection. In Fig. ^d) we enhanced 
the contrast of the these replicas for better visibility. 



We now discuss the main differences arising from the 
distortions in the surface and first few subsurface layers. 
The bandstructures corresponding to the undistorted and 
distorted 2x1 and ^/2 x ^/2 superstructures are shown 
in Fig. [Sj A striking difference between the two sets 
of calculations is that the distortions lifts the degener- 
acy between bulk and surface bands. In Fig. |9]we show 
the summary of the Fermi surfaces for various Co-doped 
slabs. The upper panel shows the result for the undis- 
torted slabs (with 2x1 and ^/2 x ^/2 Ba termination 
layers) and in the lower panel including the distortions of 
the Ba, and underlying As-Fe2-As block. The left panel 
belongs to the 2x1 and the right panel to the \/2 x 
slabs. A general result of our calculations is that the 
virtual crystal approximation at the Fe-site leads to a 
rather rigid band shift in the energy window of ±1 eV 
around the Fermi level. This is in contrast to the strong 
non-rigid band behavior upon VGA doping at the Ba-site 
[27] and can be understood as the effect of the increased 
nuclear charge at the VGA Fe- sites, which changes the 
effective electrostatic potential for the iron layers. The 
main effect of VGA Go-doping hence is an increase of the 
Fermi radii around the M-point and a decrease of the 
Fermi radii around F. This trend is seen when compar- 
ing the FS of the undoped compound in Fig. 4(b) of the 
main paper with Fig. [9]^a) andjojc). The conclusion that 
can be drawn from Fig. |9] is that the undistorted slabs 
(panel (a) and (c)) strongly resemble the /c^-projected 
bulk Fermi surfaces, except for the replica FS sheets in 
panel a, which distinguish between the x \/2 and the 
2x1 surface superstructures. Around the M-point in the 
center of the pictures the two 90° rotated ellipses, known 
from the bulk BS, are identifiable. (Note, that not all 
projected bulk bands are present in these pictures due to 
the finite slab size.) The comparison of the Fermi sur- 
faces of the undistorted (upper panels a and c) and the 
distorted (lower panels b and d) slabs shows that there is 
a strong deviation of (especially) the surface bands from 
the tetragonal symmetry due to the symmetry break- 
ing caused by the buckling and lateral distortions of the 
first few surface layers. This is more pronounced for the 
V2x V2 structure [Fig. M 



